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Abstract 

The one-dimensional Schrodinger equation for the potential + ax^ + l{l + l)/a:^ 
has many interesting properties. For certain values of the parameters I and a the 
equation is in turn supersymmetric (Witten), quasi-exactly solvable (Turbiner), and 
it also appears in Lipatov's approach to high energy QCD. In this paper we signal 
some further curious features of these theories, namely novel spectral equivalences 
with particular second- and third-order differential equations. These relationships 
are obtained via a recently-observed connection between the theories of ordinary 
differential equations and integrable models. Generalised supersymmetry transfor- 
mations acting at the quasi-exactly solvable points are also pointed out, and an 
efficient numerical procedure for the study of these and related problems is de- 
scribed. Finally we generalise slightly and then prove a conjecture due to Bessis, 
Zinn- Justin, Bender and Boettcher, concerning the reality of the spectra of certain 
T'T-symmetric quantum-mechanical systems. 
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1 Introduction 



The main subject of this paper will be the spectrum of Schrodinger equation 

r (f + 
n{a, I) il^ix) = ^-^ + x^ + ax^ + -^-^l i^i^) = E ^Pix) , (1.1) 

with boundary conditions ^(x) ^ as x — > oo and ■ip{x) ~ x'^^ as x ^ O.Q 

Over the years, many interesting features of this problem have been uncovered. For 
a < and /(/+1) = 0, the theory corresponds to a double-well potential on the full 
real line, with I = —1 and I = selecting the even and odd wavefunctions respec- 
tively. Such potentials have long been studied as toy models for instanton effects in 
quantum field theories [Q]. Furthermore, at {a = —3,1 = — 1} the ground-state energy 
£"0 is zero, with the remaining energy levels matching those at {a = 3, Z = 0}. This 
reflects the fact that for these values of the parameters the model provides one of the 
simplest examples of super symmetric quantum mechanics More peculiar are the 
properties of the spectrum at a = —(4 J + 21 + 1): a finite number of energy levels can 
be exactly computed ^ being roots of particular polynomials [^. At these points, 
the model is said to be 'quasi-exactly solvable'. Finally, at least a couple of physically 
interesting spectral problems are related to the above equation via simple variable and 
gauge transformations. The first is a linear combination of the harmonic and Coulomb 
potentials: 

2 ^ , 7(7 + 1)- 



dx"^ ~^ X ~^ 



ifix) = A ifix) , (1.2) 



with the two sets of parameters being related by 7 = {l—l/2)/2, A = —a/2, a 
2~^/'^E. The second theory is 



1 /3 5{5 + l) 



dx"^ ~^ x^l"^ ^ x^ x^ 



x(x) = Ax(x), (1.3) 



with /3 = 16aE-^, 6 = (1-3/2) /A, A = -4096^;"'^. Surprisingly, at /? = and 
6{6+l) = the latter equation is related to the Odderon problem in QCD |Q. 

In this paper we discuss some further exact spectral relationships, relating the prob- 
lems (^]^) at different values of a and /, and linking them to certain third-order differ- 
ential equations. We shall find these equivalences in the framework of the 'ODE/IM 
correspondence' [0. Some features of this correspondence are reviewed in sections 2 
and 3, and they are then used to establish five spectral equivalences in sections 4, 5, 6 
and 7. Section 8 provides an alternative insight into the second and third equivalences 
using the so-called Bender-Dunne polynomials, and in the process uncovers differential 
operators which generalise the action of the supersymmetry generators to all quasi- 
exactly solvable points. The conclusions are in section 9, and there are two relatively 
self-contained appendices. The first explains a simple but efficient numerical approach 

* The other natural boundary condition at the origin is 'tpix) ~ a;~'; when we are interested in this 
spectral problem, we shall write the differential operator as Ti.{a, —l—l). 
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to Schrodinger problems with polynomial potentials, and the second uses ideas associ- 
ated with the ODE/IM correspondence to prove a generalised version of a conjecture 
due to Bessis, Zinn-Justin, Bender and Boettcher. 



2 Bethe ansatz equations for the x^^'^ + x^^ ^ potential 

The role of functional relations in the spectral theory of the Schrodinger equation has 
been extensively explored by Voros but only recently has it been realised that 
they can lead, in certain cases, to a precise connection with the theory of integrable 
models This so-called 'ODE/IM correspondence' has been developed in a number 
of directions §-16], some of which are reviewed in In this section, we summarise 
the results obtained in this context by Suzuki concerning the Schrodinger equation 
with potential + ax^~^^ which includes the sextic potential as a special case. 



Our treatment is perhaps a marginal simplification of that given in ||T^, but the only 
genuinely new contribution to the discussion is the inclusion of the angular momentum 
term. The ODE to be considered is 



n{M,a,l)^{x) 



+ x^^' + ax''-' + ^-^l^ 



^{x) = E^{x) 



(2.1) 



with M a positive real number, which for technical reasons will sometimes be taken 
to be greater than 1. The goal is to find the eigenvalues {Ei}, those E for which 



( |2.1[) has a solution vanishing as x — > +oo, and behaving as x^ as x — > 0. (For 
3f?e/ > —1/2, the latter condition is equivalent to the demand that the usually-dominant 
x~^ behaviour near the origin should be absent; outside this region, the problem is 
best defined by analytic continuation^.) The starting-point is the uniquely-determined 
solution Y{x, E,a,l) which has the following asymptotic for large, positive x [19| : 



Y{x,E,a,l) 



-A//2-0/2 



exp 



M+1 



(2.2) 



and an associated set of functions : 
Yk{x,E,a,l) = J7'=/2+W2y(j7-fc^^ n-^^'^''E, n^^'^+^^^a, I) , n = exp(^) . (2.3) 



For integer k the Y^'s are solutions of ( |2.1| ), and any pair {Yfc,Yfc_(_i} forms a basis of 
solutions. In particular, Y^i can be written as a linear combination of Yq and Yi, the 
precise relation being 



T{E,a,l)Yo{x,E,a,l) = Y^i{x, E, a,l) + Yi{x, E,a,l) . 



The coefficient 



T{E,a,l) = W[Y^i,Yi] = Det 



y_i(x) Y,{x) 
YUix) Y{ix) 



(2.4) 



(2.5) 



^For a discussion, see chapter 4 of 
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given here as a Wronskian, is called a Stokes multiplier. In terms of the original function 
Y{x, E, a, I), ( p.4| ) taken at a and at —a leads to the following pair of equations: 

T(+)(^)y(+)(x,s) = l^-(i+°)/2y(-)(i7x,J]2ME)+l^(i+-)/2y(-)(i]-ix,f)-2^^^) (2.6) 

where 

rW(E) =r(^,±a,0, Y^'^\x,E) = Y{x,E,±a,l). (2.8) 
Keeping 3?ei > —1/2, the leading behaviour of Y near the origin at generic E is 

Y{x,E,a,l) D{E,a,l)x'^ + ... . (2.9) 

At the zeroes of D{E), the leading term is instead proportional to x'^^, and Y{x, E, a, I) 



is an eigenfunction of (2.1). This implies that D{E) is proportional to the spectral 
determinant for the problem. (For M > 1 the order of D{E) can be shown to be less 
than one, so D[E) is fixed up to a constant by the positions of its zeroes.) Setting 



D^^\E) = D{E,±a,l), from (^) and i^) we have 

T^+\E)D^+\E) = l^-(2'+l+")/2/){-)(172^^E)+l^(2m+a)/2^(~)(j^-2A^^)^ (2.10) 
T^~\E)D^-\E) = 0-(2'+l-°)/2/)(+)(f72M^)^j^(2Z+l-a)/2^(+)(5^-2M^)_ (2.11) 

Now let the zeroes of D'^^XE) be at {E^^^}, and set E = E^^'^ in either ( p^ ) or ( pl| ). 
Both T^^\E) and D^^\E) are entire in E, so the LHS of the relevant equation vanishes. 
Factorising the functions D^^^ [E) as products over their zeroes, for M > 1 the following 
system of 'Bethe ansatz' type equations for the energy spectrum is obtained: 

\ - .0-2'-!-"; (2.12) 




-0 2M^(+) j 

" , ^ = (2.13) 

Note that the spectra of the Hamiltonians TC{M,a,l) and 7i{M, —a,l) are completely 
tangled up by the Bethe ansatz constraints. 

One of the products of the ODE/IM correspondence was the realisation that energy 
levels for Schrodinger problems can be calculated using nonlinear integral equations [^. 
For the ODE {2A) with /=0 and a small, these were derived in |14|, and it is straight- 
forward to include the effect of the angular momentum term. Suppose that 

(a) all the zeroes of D^^\E) lie on the positive real axis of the complex E plane, and 

(b) all the zeroes oiT^±){E) lie away from it. 

(As shown in appendix B, (a) holds if / > -1/2, and (b) if |a| < M + 1 - |2/+1|.) 
Setting 

^ ' L>(=F)(i^2A/^) ^ \M+i y ' ^ ' 
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the equations are then 



i-{2l + 1 



a , 
M' 



ib, 



'0 e 



+ 



de'Ki{e 
dO' K2(e 



3'- 



ln(l + a±(0')) - I d9'Ki 

JC2 



- e') ln(l + 

e')\n{l + a^{e'))- / dO' K2{e - e')\n{l + 

C2 



{6') 



-) (2.15) 



The integration contours Ci and C2 run just below and just above the real axis respec- 
tively, and the constant 60 = 7r-'^/^r(2^)/(2Mr(| + )) is fixed via a consideration of 
the WKB asymptotics of D^^\E) for \E\ — > 00, arg(£') 7^ 0. An integral expression for 
the kernels Ki and K2 at general values of M was found in |14]. In this paper we are 
particularly interested in the sextic potential, and here we note that the special nature 
of this case is reflected in the fact that at M = 3 the kernel functions can be explicitly 
integrated, and have the simple forms 



V3 



K2{e) 



27r(2cosh6' - 1) 



(2.16) 



27r(2 cosh (9 + 1) 

A search along the real axis for the zeroes of the functions 1 + a^{0) provides the 
energy levels of the Hamiltonians Tl{M, ±0,1). While we do not have a rigorous proof, 
we expect that the solution to ( |2.15| ), which can readily be obtained numerically by 
iteration, is unique for a and / in the stated range. This is one way to justify the claim 
that, for / > —1/2 and \a\ <M + 1 — |2/— 1|, the Bethe Ansatz equations ( 2.12| ) and 
(|2.13D , together with the 'analytic properties' (a) and (b) and the WKB asymptotic 
which determined 60, characterise the set of numbers {E^^} uniquely. (An alternative 
approach to this question might be to generalise the analysis of @,H and appendix A 



of |20|, based on the so-called quantum Wronskian relations.) 

To treat more general values of a and I, we will find it most convenient to work 
directly with the Bethe ansatz equations. In contrast to the integral equation, these do 
not have a unique solution. As is standard in studies of the Bethe ansatz, it is useful 
to take logarithms. For / > —1/2 and \a\ < M + 21 + 2 (see appendix^, all energies 

are positive. For this situation we assume that there is a 



E 



(±) 



E, 



(±) p{±) ^(±) 







^1 



one-to-one correspondence between these energies and the integers /c = 0, 1, 2, 




-ZTT 



-ITT 



21 + 1 + a 
M + 1 

21 + 1- a 



+ 2A; + 1 



M + l 



+ 2k + l 



(2.17) 



(2.18) 



where the logarithms are all on the principal branch: — vr < — i In < vr. At larger values 
of \a\ and \l\ some of the low-lying energies might become negative, and in such cases 
care must be taken to keep track of the nontrivial monodromy of the log function. 



4 



3 The Bethe ansatz approach to a third-order equation 



This section summarises the derivation of equations of Bethe ansatz type for a third- 
order ODE with 'potential' x^^, fohowing |12,TH]. We start with the equation 



ax'^ x"^ \x^ ax x^ / i 



L 

X 



1 d 1 

. x"^ dx x^ 



and, as in |15|, rewrite it as 

[D(g) + x3^](/<(x) = i?(/>(x) 
where g = {go-, 91,92} with 90 + 9i + 92 = 3, and 

D(g) = D(g2 - 2)D(5i - l)D(5o) , D(g) = 



d 9 



(3.1) 



(3.2) 



(3.3) 



, dx x , 

The relationship between g and {G, L} is 

G = 9o9i + 9092 + 9192 - , L = 2 - 909192 - {9091 + 9o92 + 9192) ■ (3.4) 

Again we introduce a uniquely-defined function y{x,E,g), which solves ( |3.1| ) and tends 
to zero as X ^ cxD along the positive real axis as 



-N 



y{x,E,g) 



,JV+1 



expl-^ 



(3.5) 



Given y{x), bases of solutions are constructed just as in the second-order case. Set 



yk{x,E,g) =uj''y{uj ^x,uj ^'"^ g) , w = exp 



3Nk 



2m 
3Af+3 



(3.6) 



For integer k, yu solves (0.1D , and {y^, yk+i, yk+2} form a basis. We can therefore expand 
yo as 

yo - C^'hi + C(2)y2 - ys = (3.7) 

with coefficients C^^^ and C^^-* - Stokes multipliers - which are independent of x. Elim- 
inating C^'^\ we have 

W02 - C(i%i2 + W23 = , (3.8) 
where the Wronskians of pairs of solutions, 



Wkik2 = W[yk^ yk^] = Det 



ykiix) ykiix) 
y'kS^) y'k^i^) 



(3.9) 



were used. (Note, since y^-^ and solve a third-ordei equation, these are nontrivial 
functions of x.) Now multiply by yi and use the relation yiWo2 = yoWi2 + 2/2l^oi to 
find 

C^'hiWu = yoWu + y2Woi + yiW23 . (3.10) 
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In this form the relation can be rewritten in terms of just two functions, y{x, E, g) and 
Wix,E,g) = Woi{x,E,g) : 

C'^^\E)y{uj''^x,uj~''^^ E)W{uJ-'^x,uj-^^E) = uj~^y{x,E)W{uj'^x,uj-^^E) + 

y{uj-^x, uj-^^E)W{x, E) + ujy{oj-^x, uj~^^ E)W {uj'^ x , E) . (3.11) 

We initiahy suppose that 5Re((7o) < '^^{gi) < 3^e((72)- Then as x — > the leading 
behaviours of y and W are 

y{x) ~ D^^^ {E, g) x^o , W{x) ~ Z)(2) (E, g) x^o+^i"^ . (3.12) 

Using equations ( p.l2| ) the relation ( p. 11] ) becomes 

C(^)(i?)Z)(l)(c^-3^i^)Z)(2)(c^-3JV^) ^^90-1^(1) (^)^{2)(^-37V^) ^ 

^9i-lZ)(l)(^-6^i?)Z)(2)(i?) +,,2-,o-9iZ){l)(^-3iV^)^(2)(^-6iV^)_ (3.13) 

Again we shall consider this expression at the zeroes of D^^^ and D^'^\ It is convenient 
to write these as E^^^^ and uj'^^^'^E^jf' , so that 

Z)a)(i?«,g) = 0, D^'){u:^^l'Ef\g) = 0. (3.14) 

As in the second-order case, these functions have a spectral interpretation. In particular, 
the vanishing of D^^\E) signals the existence of a solution, at that value of E^ to ( |3.2| ), 
decaying as x ^ oo, and having a faster-than-usual decay at the origin: 

y{x) ^ x'^'''^3i,92) ^ 3.^0. (3.15) 

(Since gi and g2 can be complex, by min((/i, (72) we mean whichever of gi and g2 has the 
smallest real part.) Evaluating ( PH) at ^ e {^^i^^} and E £ {lo'-^^/^ Ej^^} and imposing 
the entirety of C^^\E) leads to the following set of S'C/(3)-related BA equations, with 
A; = 0,1,2,...: 




eP -u-^e^^ ' 

- E^P 



-Lo3°-a^ ; (3.16) 



^1 = -,^291+50-3, (317) 



Since (70+51+92 = 3, the right-hand sides of equations (|3.16| ) and ( |3.17|) can be given 
a more symmetrical appearance by rewriting them as — ct;2go+S2-3 g^j^^ _^;-292-3o+3 
respectively. These equations, together with WKB-like asymptotics for D^^\E) and 
D^^\E), fix the numbers E^^ and E^^ up to discrete ambiguities, which for the prob- 
lems in hand can be eliminated given some facts about the approximate positions of 
the zeroes of the functions D^^''^\E) and some associated functions a^^''^\E). These 
are analogous to the analyticity conditions (a) and (b) of the previous section, and are 
described in more detail in |T2| , p!5[| . It is also possible to solve the system via a nonlinear 
integral equation, but this will not be needed here. 
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4 The first spectral equivalence 



The first spectral equivalence follows from observing that at = 1, a; 

3N 

io~ = i. The SU{3) BA equations therefore simplify to 



3N 




-ijj 



-ijj 



290+32-3 . 



-232-90+3 



-1 and 



(4.1) 



(4.2) 



These equations coincide with the system ( p. 121) , ( p. 13 ) at M = 3 provided the right- 
hand sides of the two BA sets are equated: 



(250 + 52 -3)/3 = (-2/- l-a)/4, - (2^2 + 5o - 3)/3 = (-2/ - 1 + a)/4 . (4.3) 

Combined with a matching of the analytic properties (a) and (b), this suggests the 
following relationship between quantities in the two problems: 

D^'HK-'E,g) = f{a,l)D^+\E,a,l), 

D^^'>{iK~^E,g) = f{a,l)D^-\E,a,l). 



(4.4) 
(4.5) 



The proportionality factors f{a, I) and k cannot be determined by a comparison of the 
Bethe ansatz equations alone. However, as in [12|, k can be calculated by comparing 
the large negative E asymptotics of D^^^ and D^. The result, independent of a and I, 
is K = 4/(3\/3). Solving (|4.3| ), the parameters {a, I) and g are related as 



a 



0(50,52) = 2(2-50 -92), 1 = Kgo,92) = (252 



25o)/6. 



and 



50 = (l-a-6/)/4, gi = (l + a/2) , g2 = {7 - a + Ql)/4 . 
Thus we have a spectral equivalence between the following eigenvalue problems 



6 2 + 
+ x^ + ax^ + ^ ^ ' 



K 

where 



dx'^ 
G = 



x"^ dx 



<I>(x) 

(t>{x) = 



E 



E $(x) , 
x) , 



X 



mm(9i,92) 



(4.6) 
(4.7) 

(4.8) 
(4.9) 

(4.10) 



5o5i + 5052 + 5152 - 2 , L = 2 - 505152 - (5o5i + 5o52 + 5152) , 

and the parameters in the two models are related as in (1^ and (13). Note that the 
general 5C/(3)-related equation at = 1 is mapped onto the general sextic-potential 
problem. The number of parameters matches up, because the third-order equation 
allows for two linearly-independent angular momentum type terms. The different roles 
that these parameters play in the second-order problem will be important for the next 
spectral equivalence that we discuss. 
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5 The second spectral equivalence 



Our second spectral equivalence is related to the enhanced symmetries of the third-order 
equation. To be more precise, the differential equation (3.2) is unchanged under per- 
mutations of {90,91,92}, while the values of a and / which appear in the corresponding 



second-order equation, given by (4.6), are not. If we make a continuation in {50)51)52} 
which swaps 51 and 52 while leaving 50 unchanged, then both the third-order equation 
itself, and the specification ( 3.12D of D^^\ are unchanged, and so 



D«(K-ii?,g)^Z)«(K-ii?,g). 



Using (4^), this means that 



f{a,l) 



(5.1) 



(5.2) 



with 



5 = a(5o,5i) = (3-a + 60/2 , / = /(50, 5i) = (a + 2Z - l)/4 . (5.3) 
It will sometimes be convenient to put this in matrix form. If a = {a, 1, 1)"^, then 



-1/2 3 3/2 
a = Ta , with T = ( 1/4 1/2 -1/4 

1 



(5.4) 



At this stage we do not know how to calculate /(a, /) and /(a, /) exactly, but (|5.2| ) can 
be combined with (B.15) from appendix B to give their ratio: 



/(a,/) _ r(^ + i) 



f{a,l\ 



(5.5) 



Note that this is singular or zero at negative-half-integer values of / or at which a 
'resonance' is expected in one or other of the spectral problems ||T^. Away from these 
points, we have a spectral equivalence between 



6 2 + 
+ x° + ax^ + ^ — ^ 



^x) = E^{x) , «>U^o~x'+i 



and 



dx"^ 



(3-a + 60 2 {a + 2l-l){a + 2l + ?,) 



-x^ + 



16x2 



(5.6) 



(5.7) 



with <^|a;^o ~ 2;(a+2i-i)/4+i^ alternative viewpoint on this equivalence in terms of 
intertwining operators will be given in §8 below, while some direct numerical checks are 
reported in appendix A. 



8 



6 The third spectral equivalence 



As was mentioned in the introduction, at the special values a = aj{l) = — (4J + 2/ + 1), 



with J a positive integer, the model (1.1) is 'quasi-exactly solvable' (QES), and the 
first J energy levels can be computed exactly. For J = 1, the single exactly-solvable 
energy is the ground state and the model is an example of supersymmetric quantum 
mechanics. This is signalled by the fact that the Hamiltonian at a = ai{l) = — (2/+5) 
can factorised in terms of first-order operators as 



H{ai,l) 



(6.1) 



where 



Q- 



ax 



l + l 



-r + x-" 
ax 



l + l 



(6.2) 



The SUSY 'partner' Hamiltonian 7i = 7i{ai,li) is obtained through the intertwining 
relation W(ai,Ti) Q+ = Q+ n{ai,l) with ai{l) = 1-2/ and Ti{l) = l+l : 



n = Q+Q- 



(6.3) 



The wavefunctions of the two models are simply related by 



V'j(a;) = Q^^i(x) 



(6.4) 



However, the ground-state wavefunction of H{a,l) is -00 = x'+-^exp(— ^) , and this is 
annihilated by Q"*". As a result, 7i and 7i are spectrally equivalent save for the extra 
level at E' = only present in Spec('H). This (very standard) result makes it natural 
to ask whether similar 'partner potentials' might exist at higher values of J, sharing 
the same spectra as the QES Hamiltonians 7i{aj{l),l) apart from the first J levels. We 
shall find that this question has a surprisingly simple answer by using the Bethe ansatz 
approach to the spectral problem. 



Setting a = aj{l) = — (4J+2/+1), the J exactly-solvable levels Eq~^\ e[^\ 



...E 



lie in the sector '(+)'. The BAE ( ^TzD , (|2T^ ) for M=3 are then 



{+) 
J-i 



J-i 

n=0 



>(+) 



lE, 



(-)' 



Ei+'^ + iE, 



(-) 



n=0 

oo 

+ ^ln 

n=J 



E: 



(-) 



iE 



E.. 



(+) 



iE 



-i7r[-J + 2k + l] 



(6.5) 



-i7r[l + 3/2 + J + 2k], (6.6) 



where the integer k runs from to cxd. Next, we will use the fact that the exactly- 



solvable energy levels appear symmetrically, as E^ 



(+) 



-E 



(+) 



-1' 



to simplify the first 



sum on the LHS of 



Recalling that, since a is negative, the E^ ^ are positive for 
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/ > —1/2, and keeping track of the monodromy of the logarithms by using the reflection 
formula 



In 



we obtain 



-X + i 
J-i 



-27ri — In 



n=0 



X — I 

X + i 
(-)\ 



(x > 0) , 



-ivr J . 



(6.7) 
(6.8) 



Finally, ignoring the first k = 0, 1, . . . J— 1 instances of ( ^.5| ) and relabelling E^j^j 
E^^^ we end up with 




-m[J + 2k + l] ; 
+ 2k + 3/2] 



(6.9) 
(6.10) 



where k again runs from to oo. Comparing with ( p^ ), ( ^Isj) we now reinterpret 
the left-hand sides of equations ( |6.9D and ( |6.10| ) as the quantisation conditions for the 
energy levels of a new potential, with parameters aj and Ij : 



- z7r[J + 2A: + 1] 
-m[l + 2k + 3/2] 



-ZTT 



2/j+l+Sj 



2l,j+l-aj 



+ 2A: + 1 



+ 2A: + 1 



Solving, 

Sj = 2J- 2/ - 1 = -(3 + a + 60/2 , ^J = J + I = {-a + 21 - I)/ A, 
and this can again be put in matrix form, for a = (a, 1, 1)^, as 



a 



. a 



with 



-1/2 -3 -3/2 
4/4 1/2 -1/4 
1 



Then, for J G N, 

= 4^^(aj,t/) , (A: = 0, 1, 2, . . .) 
and, modulo the exactly-solvable levels, there is a spectral equivalence between 



(6.11) 
(6.12) 

(6.13) 
(6.14) 
(6.15) 



+ - (4J+2m)x^ + ^ 



dx"^ 



$(x) = £;$(x) , 



(6.16) 
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and 



a , , ^ 2 {J+l){J+l + l) 
--^+x^ + {2J-2l-l)x^ + ^ 



^{x) = E , $U_.o ~ . 

(6.17) 

If J is a negative integer, the mapping still makes sense, but it acts in the opposite sense: 
shifting / — > /—J, (|6.17 ) becomes the QES problem for a\j\^ and ( 6.16 ) its partner with 
the QES levels removed. 

Finally, we still have the freedom to apply the 'tilde-duality' T discussed in §5, so 
( |6.17D is in turn isospectral to 



x^ 



^x)=E^{x), $|a;^o~a; 



J-l/2+1 



(6.18) 

This chain of equivalences will be discussed further in the conclusions. 

Since the quasi-exactly solvable energies and the associated wavefunctions are in 
principle exactly known, one could eliminate them one by one using the Darboux trans- 
formation, though this would be a lengthy business for large values of J. What seems 
surprising about the results ( |6.17 ) and ( 6.1§| ) is that the potential can have such a 
simple form once all of these levels have been subtracted. 



7 The fourth and fifth spectral equivalences 

The equivalence of the second-order equation with an 5'C/(3)-related third-order equa- 
tion suggests two further spectral equivalences. As explained in []T2| , [T5| , the Z2 symmetry 
of the SU (3) Dynkin diagram is reflected in a relation between the functions y and W 
that were introduced in §3. Explicitly, 

y(x,ii;,gt) =Ty[y_i/2,yi/2](a:,i?,g), (7.1) 

where = {(^q, g|, 52) and g\ = 2—g2_^ ■ On the second-order side of the story, a similar 
Wronskian appears, but this time in the formula (|2.5D for T{E,a,l): 

T{E,a,l) =W[Y^uYi]{E,a,l). (7.2) 

Before the two equations can be compared, the x-dependence must be eliminated from 
( [7.ID , and as usual this is done by considering the behaviour as x — > 0. Extending the 
definition (|3.12| ), we define three functions ^ with D^^^ = D^q^I , by 

2 

y{x,E,g)=Y^ ^[^ (E, g) Xi {x, E, g) , (7.3) 

where the solutions Xiix,E,g) to ( |3.2| ) are defined by Xi ~ x^' + 0(a;^''^^) as x ^ 0. 
To match ( |7.2[ ), we expand out the RHS of ( f7.1| ) and then project onto the component 
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behaving as x^"^^'^ = x^i to find 



On the other hand, the function Y{x, E, a, I) can be expanded as 

Y{x, E, a, I) = D{E, a, l)X{x, E, a, I) + D{E, a, -l-l)X{x, E, a, -l-l) 



(7.4) 



(7.5) 



with X{x, E, a, I) ~ x~^ as a; ^ 0. (Cf. eq. (5.2) of [T^, but note that the definition of 
D{E, I) used in differs from that used here by a factor of {2l+l)~^. ) Now substitute 
into (7.2) taken at {—E, —a, I) : 



T{-E,-a,l) = {21+1] 



-21-1 



D{-n'^'E, a, l)D{-n-'^' E, a, -l-l) 



21+1 : 



D{-n-^^^E,a,l)D{-n^"'E,a,-l-l) . (7.6) 



2M ■ 



2M 



UJ' 



3N/2 _ g7i-«/2^ and, if g and / are related by 



If M = 3 and TV = 1, then -Q 

^{5o-92)/2 ^ ^-21-1 ^ g-7ri(2«+i)/4^ Furthermore, from (U), 

D[^J{E,g)=f{a,l)D{KE,a,l) , D[lJ{E,g) = f{a,-l-l) D{KE,a,-l-l) . (7.7) 

(The second relation is obtained by a continuation which swaps go and g2-) Using these 
identifications and comparing (7.4) and ([7^) gives our fourth spectral equivalence: 



Z)f.V(K-ii^,gt) = :^f{aJ)f{a,-l-l)T{-E,-a,l). 



[1] 



(7.8) 



where gj = 2 — g2_^ and a = 2(2 — go — g2), I = (252 — 3 — 2go)/6. As with the 
first equivalence, this relates spectral data for differential equations of different orders. 
The spectral interpretation of functions such as T in the ODE/IM correspondence was 
discussed in |11|, and is reviewed and extended to the current context in appendix B 
below. 

We can obtain a relation between objects in the second-order equation by using the 
first spectral equivalence to rewrite the LHS of ([7.8D, and this constitutes our fifth and 



final spectral equivalence. The only subtlety is that (}4.4D involves d'^q'I , not d'^I'I , and 
this can be overcome by a continuation in the gi . Swapping go and gi and tracing back, 

2/(-5,r 



T{-E, -a,l) 



3/(a, /)/(«,-!-/) 



D{E,-a{a,l),l{a,l)) . 



(7.9) 



The proportionality factor can also be found explicitly, by considering (7.9) at E = 
and using formulae (B.15) and (B.16). The result: 

2 \/f7r ~ 
T{-E,-a,l) = —D{E,-a{a,l),l{a,l)). (7.10) 



r(/(a,o + 
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Via the second equivalence, this can be rewritten as 



2 \ m ~ 

T{-E,a,l) = ^ — D{E,a{a,l),-l-l{a,l)) . (7.11) 

T{-l{a,l) - 2 ) 

As explained in appendix B, T is the spectral determinant for a 'lateral connection' 
problem, with the wavefunction lying on a contour in the complex plane joining a pair 
of Stokes sectors at infinity. In contrast, D is the spectral determinant for a 'radial 
connection' problem, with the wavefunction living on a half-line. Similar equivalences. 



albeit for slightly different potentials, have been found in |21|, and it would be interesting 
to see whether similar methods could be applied in this case. 

The mappings of parameters involved in these relations can be streamlined by in- 
troducing two further matrices, A and L, again acting on the vectors a = {a, 1, 1) 



T. 





r-1 





0^ 









1 









^ 





1 ) 







-1 -1 I , (7.12) 


and setting 

7(a) = 7((a, I, if) = 2^/r(/+i) . (7.13) 



The tilde-duality (5^) is then 

-f{a)D{E, a) = -f{Ta)D{E, Ta) , (7.14) 

while (^J^) and ([7Tl|) are, respectively, 

T{-E, Aa) = -/{ATa)D{E, ATa) , T{-E, a) = -/(LTa) D{E,]LTa) . (7.15) 

Note also that EI = ATA, so the first relation of ( |7.15| ) can be rewritten as T{—E, a) = 
^ (Ma) D{E, Ma). At the QES points, a = aj = {aj{l),l,l)'^ , the spectrum encoded 
by D{E,Ilaj) is equal to that of D{E,aj), apart from the QES levels. As will be 
explained in the next section, these levels are in fact the zeroes of the Bender-Dunne 
polynomial Pj{E), so at the QES points we have 

Pj{E)T{-E, aj) oc D{E, aj) , (7.16) 

which is a relation between spectral data for lateral and radial connection problems 
with the same Hamiltonian. 

The algebra of the matrices we have introduced is best described by first defining 
M = AL. Then a set of defining relations for L, M and T is 

= = = (LM)2 = (MT)2 = (LT)^ = I . (7.17) 

Thus M commutes with L and T, while < L, T > forms the Weyl group of SU{3). 
However, in general only T yields a spectral equivalence of the D{E, a). We will return 
to this point in the conclusions. 
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We end this section with two further remarks about the fifth set of equivalences. 
First, they can also be obtained entirely in the context of the second-order differential 



equation. For M = 3, manipulating equation (2.11) and using the fact that Q 



)2M 



leads to the following functional relation, special to this particular value of M: 
2sin(f (2/+l-a))L'(+)(^) = 

j^(2Z+l-a)/2^(-) .^)^(-) (_ -^^ _ ^-(2i+l-a)/2^(-) ^iE)D^-^ (iE) . (7.18) 



Taking (|7.18|) at ^ = E^^^ , combining it with (|2.1l|) , also at ^ = eI~^\ and finally 



.(+) 



expressing the result in a factorised form over the zeroes of T(-)(£;) (which we denote 
by {—'^k '*} ) yidds the following set of constraints: 



n 

n=0 



An 



A complementary set arises from ( p. 10 ), taken at E 



n 

n=0 



-C/n 



-2AI ^{-) 



ii;(+)_0 2M^(-) 



-n 



2i+l-a 



(7.19) 



(7.20) 



Together, ( |7.19D and ( |7.20D form a set of Bethe ansatz equations of exactly the same 
form as ( |2.12|) an d ( 2.13| ), save for the replacement of a and / on the right-hand sides of 
( |2.12D and ( |2.1^ ) by a{a,l) and l{a,l), respectively. By comparing the left-hand sides 
and exploiting the analytic properties derived in appendix B, one obtains by another 
route the second and fifth equivalences (^), ([7.9^: 



E^^\a,l)=E):>{a,l) and Xr\a,l) = E)~' {a,l) . 



.{+) 



(7.21) 



The previous approach, which proceeded via the symmetries of the third-order equation, 
was perhaps more elegant. The advantage of this alternative method is that the only 
analytic properties used are those of spectral determinants of the second-order equation, 
and these are known rigorously from the results in appendix B. 

The second remark relates to the fact that the lateral connection problem solved 
by T is closely related to PT-symmetric quantum mechanics. These problems are not, 
in any obvious sense, self-adjoint, and the reality properties of their spectra have long 
been of interest |22-^. Since the spectral problems on the right-hand sides of ( [7.10 ) 
or (7.11) are self-adjoint for I > —1/2 (respectively —l—l > —1/2), these two identities 



give us a simple way to understand the reality of the spectra encoded by T in these 
particular cases. However, in appendix B below we will give a proof of reality which is 
both more general and more direct, and so we will not pursue this any further. 
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8 A relation with Bender-Dunne polynomials 



The dualities that we have been discussing have an interesting relationship with the 
so-called Bender-Dunne polynomials. These were introduced in Q as a way of under- 
standing quasi-exact solvability, but here we will also be interested in their properties 
at general values of the parameters. Briefly, one searches for a solution to ( |1.1| ) of the 
form 

*,.,E,a.0^e--/^.«f (8.1) 

n=0 

For this to solve the differential equation ( |1.1D , the coefficients P„ must satisfy the 
following recursion relation: 

Pn{E) = EPn^i{E) + 16{n-l){n-J-l){n + l-l/2)Pn-2{E), (n > 1) (8.2) 

where, as before, J = J{a,l) = — {a+2l+l) / 4. The value of Po{E), which determines 
the normalisation of ip{x)^ is conventionally taken to be 1; from (|8.2| ), Pi = E, and P„ 
is a polynomial of degree n in E, known as a Bender-Dunne polynomial. So long as 
/ 7^ — n— 3/2 for any n G Z^, ( ^.1| ) will yield an everywhere-convergent series solution 
to ( ll.ll ). Furthermore, this solution automatically satisfies the boundary condition 
ip ~ x'"*"^ at X = ; but at general values of E, it will grow exponentially as x — > oo. We 
now ask whether there are transformations of the parameters a and / which leave the 
Bender-Dunne polynomials invariant. It is easily seen that if J and / are replaced by 
J = —1 — 1/2 and / = — J — 1/2, then the recursion relation is unchanged. Translated 
back to the parameters a and I, this implies that 

PniE,a,l)=Pn{E,aJ) (8.3) 

where 

5 = 2J + 4/ + 2 = 3/2-a/2 + 3/ , 7= - J - 1/2 = q/4 + //2 - 1/4 . (8.4) 

This matches the 'second spectral equivalence' found earlier. 

We will return to this case later, but first we discuss the special points where the 
model is quasi-exactly solvable, for which a similar game can be played. If a and I are 
such that J{a, I) is a positive integer, the second term on the RHS of (|8.2D vanishes at 
n = J+1, and all the subsequent P„ therefore factorise: 

Pn+j{E, aj, I) = Pj{E, aj,l)QniE, aj, I) , (n > 0, J = -(aj+2/+l)/4 G N). (8.5) 

Hence, if Pj{E) vanishes then so do all Pn>j{E) and the series ( ^.1[ ) terminates, au- 
tomatically yielding a normalisable solution to (|l.l| ). The J zeroes of Pj{E) are the 
J exactly-solvable levels for the model and, as observed by Bender and Dunne, this 
provides a simple way to understand the quasi-exact solvability of the model. Now we 
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would like to go further and discuss the remaining levels. The polynomials Qn satisfy 
the recursion 

Qn{E) = EQn-i{E) + 16(n + J - l)(n - l)(n + J + 1- l/2)Qn-2{E) , (n > 1) (8.6) 

with initial conditions Qo = 1, Qi = E. This matches the recursion relation for Pn{E), 
so long as J and / in (8.2) are replaced by J = — J and I = J + 1. Hence, if 

Qj = 2J-2/-l = -aj/2-3/-3/2 , Tj = J + / = -qj/4 + //2 - 1/4 , (8.7) 

then 

Qn{E, aj, I) = Pn{E, ajjj) . (8.8) 

This corresponds to the 'third spectral equivalence', and it has an interesting conse- 
quence for the series expansion (|8.l| ), which we rewrite as 



ip{x,E,aj,l) 



_|_ ^ (-i)" -^"(-^i ^2n 



n=J 



n\T{n+l+3/2) 

/ i.n+j Pj{E,a,j,l)Qn{E,aj,l) 2 
^' (n+J)!r(n+J+/+3/2) 



(n+J) 



n=0 



(8.9) 



the dots standing for lower-order terms. This can be compared with the expansion of 
the wavefunction ij){x,E,aj,lj). Using Ij = J+l and the equality 



this is 



ip{x,E,aj,lj) 



n=0 



Qn{E,aj,l) 
n!r(n+J+/+3/2) 



„2n 



X 



^.10) 



It is now easy to see that ^{x,E,aj,l) is mapped onto a function proportional to 
ip{x, E,aj, Ij) by the differential operator 



g-a;4/4^/+J+l 



1 d 

X dx 



1 d 

X dx 



+ x' 



l+l 
x^ 



^.11) 



This is enough to see that the following intertwining relation between differential oper- 
ators must hold: 

Qj(l) Hiajil), I) = H{aj{l)Xm Qj{l) ■ (8-12) 

(Consider|| the difference between the LHS and RHS. This is a linear ( J+2)*^-order 
differential operator, independent of E, and it is easily seen that it annihilates the func- 
tions 'ip{x, E,aj{l),l). These functions are linearly independent for different values of 
E, while a (J+2)*'^-order operator can annihilate at most ( J+2) independent functions, 
unless it is identically zero. This establishes the equality.) We also used Maple to 
verify ( ^.12 ) directly. Finally, Qj{l) respects the boundary conditions: if ^(x) decays 

*We would like to thank Peter Bowcock for a discussion of this point. 
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as x — > oo then so does Qj{1)%Ij{x) , and if i^ix), given as a series by (|8.1|), has lead- 
ing behaviour x''^^ at the origin, then Qj{l) ipi^) has leading behaviour Thus 
Qj{l) maps eigenfunctions of the problem Ti.{aj,l) to those of Ti.{aj{l),lj{l)), or to 
zero. The eigenfunctions mapped to zero are those for which Pj{E,aj,l) vanishes (the 
lower-order terms '. . . ' in ( |8.1D are clearly annihilated by Qj(Z)), and these are precisely 
the exactly-solvable levels. This provides an alternative derivation of the duality found 
with the aid of the Bethe ansatz equations in §6, and shows that in Qj we have found 
the generalisation of the supersymmetry operator Q"*" = to the QES problems ( |6.16 ) 
with J > 1. 

So far we have discussed the action of Qj{l) on solutions to the spectral problem 
with x'^^ boundary conditions. But given the intertwining relation ( |8.12 ) it is natural 



to look for an action on solutions satisfying the other, x~ , boundary condition at the 
origin. As a relation between differential operators, the fact that 7^(0, b) = Ti^a, —1—6) 
means, trivially, that 

Qj{l) n{aj{l), -1-0 = n{aj{l), -1-Tj{l)) Qj{l) . (8.13) 

It can also be checked that, in general, the relevant boundary conditions are respected, so 
that (|1|) holds as an intertwining relation between eigenvalue problems. Substituting 
— l—l for / throughout, the relation is 

Qj{-l-l) n{aj{-l-l),l) = H{aj{-l-l), Qj{-l-l) . (8.14) 

Thus Qj(— 1— /) and its adjoint intertwine between the spectral problems 

W(-4J+2/+l,/) and W(2J+2/+l, /- J) , (8.15) 

and, in general, no eigenfunctions are annihilated. Furthermore, the mapping ( ^.15]) 
is exactly the 'second spectral equivalence' of §5 above, specialised to cases where the 
initial pair of parameters (a, I) satisfies a = — 4J + 2/ + 1. This hints at an alternative 
way to obtain ( ^.14|) : just as was done at the QES points using ( ^.81 ), one can compare 



the series expansions using (8.3). At a formal level, for any value of a and / the series 



for ip{x, E, a, I) is mapped onto that for ^(x, E, a, I) by 

VJl) = e--'/^x^-' ( -4-]' e^'/'x' , (8.16) 
\x ax J 

where ^ = —a/4 + 1/2 + 1/4, and the action of a fractional power of the derivative is, 
again formally, defined by 

1 dY ^2n _ r(n+l) 2(n-/.) 



x"" = 2^ ^/ ' \ x"^"-^^ (8.17) 
X ax J r(n+l— /i) 

In cases when /i is a positive integer, V^{1) becomes an ordinary differential operator, 
equal to Q^{—l—l), and ( 8.14| ) is recovered. It is interesting to speculate about the 
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existence of some kind of spontaneously-broken 'fractional' generalised super symmetry 
lying behind the tilde-duality at arbitrary values of the parameters, but we leave this 
for future work. 

In concluding this section, we would like to mention the recent article |33|, which we 
noticed as we were finishing the writing of this paper. By a completely different route, 
involving a study of a concept called 'AA-fold supersymmetry' |U, 35 1| the authors of 
this work have also introduced higher-order analogues of the supersymmetry generators. 
Although the connection with quasi-exact solvability is not mentioned, one can check 
that the 'cubic' case of the type A AA-fold supersymmetry of ||3^ reproduces the result 
( |8.12D above, albeit with a slightly different presentation of the operators. (In fact, using 
the most general form of their operators, one can also obtain an intertwining relation 
for the more general QES sextic potentials involving an additional x'^ term.) We should 
also mention that a connection between certain other forms of non-linear supersymmetry 
and quasi-exact solvability has recently been pointed out in |^^. However the forms of 
the supersymmetry generators explicitly treated in that paper do not cover the case of 
the sextic potential discussed above. 



9 Conclusions 

The main purpose of this paper has been to illustrate how spectral properties and 
symmetries of interesting differential operators can be handled using tools originally 
developed in the context of integrable models. These lead to some novel spectral equiv- 
alences, and, as will be shown in appendix B below, they also allow for an elementary 
proof of a reality property which has been surprisingly elusive when studied by more 
conventional methods. Some of the equivalences we have subsequently been able to 
re-derive by other means, and in this respect, the role of higher-order generalisations of 
the supersymmetry operators is particularly intriguing, especially in the light of their 
independent appearance in |3^. It would be very interesting to find out whether the 
connection between such operators and quasi-exact solvability that we have observed is 
more general. In this paper, we obtained the operators Qj through a direct examination 
of power series solutions; how they fit into the more algebraic schemes for understand- 
ing quasi-exact solvability, as developed in, for example, |^,^, is another question that 
deserves further study. 

In many ways, the first and fourth spectral equivalences, between second- and third- 
order equations, are the most unexpected of our results. They can be traced back to the 
collapse of the SU{3) Bethe ansatz equations at A^=l. A similar phenomenon occurs 
in 5C/(n)-related BA systems for n > 3, which are related to higher-order differential 



equations via the ODE/IM correspondence 13, 1^. However, in these cases the 
resulting 'reduced' systems are not so readily identifiable, and so at this stage we lack 
an interpretation of the phenomenon in terms of the properties of differential equations. 



^generalised supersymmetries of this sort are also called 'higher-derivative', or 'non-linear' 
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On the other hand, via the first equivalence we do at least see that quasi-exact solvability 
is not restricted to second-order spectral problems. 

One can also ask whether quasi-exact solvability might have a role on the 'integrable 
models' side of the ODE/IM correspondence. Bethe roots correspond to zeroes of D{E), 
and at the QES points the locations of a finite subset of these can be found exactly. 
Prom the identity ( 7.16| ), the locations of the remaining roots coincide with the zeroes 



of T{—E), while, as follows from the form of the Bender-Dunne wavefunctions, the QES 
roots themselves correspond to coincidences in the locations of zeroes of D^^\E) and 
r(=F)(S). However, we do not know of any special significance of these facts. 

Finally, it is interesting to draw the full set of spectral problems that can be reached 
from a QES starting-point using the dualities that we have been discussing: 

n{-4:J-2l-l,l) — > n{2J-2l-l,J+l) 



n{2J+Al+2,-J-\) — > n{2J+Al+2,J-\) 

Vertical arrows correspond to the second spectral equivalence T, while the upper 
horizontal arrow is the level-eliminating third equivalence, H. The two problems on 
the bottom row are related by the transformation L: / They correspond to 

the same Schrodinger equation, and differ only in the boundary condition imposed at 
the origin. It follows from the diagram that the 'regular' eigenvalue problem for this 
equation, that with the x"^"^/^ boundary condition, has exactly the same spectrum as 
the irregular problem, with the exception of the first J levels. This can be understood 
by noticing that, at / = J— 1/2, there is a 'resonance' between the regular and irregular 
solutions of the Schrodinger equation (see, for example, |]l8|). The Bender-Dunne series 
expansions for the first J levels of the irregular problem truncate before the resonance is 
reached, whilst the wavefunctions for the remaining levels are completely dominated by 
the effect of the resonance, and hence match those of the regular problem. In fact, such 
a square of spectral problems can be drawn starting from any values of the parameters 
a and Z, on account of the identity EI = TLT. But it is only at the QES points that the 
horizontal directions correspond to (partial) spectral equivalences, since the resonance 
between regular and irregular solutions just described only occurs when J is an integer. 
Thus we have some novel points at which the sextic potential can be considered to be 
quasi-exactly solvable, a 'dual' interpretation of quasi-exact solubility for this model 
in terms of the resonance of irregular solutions of the Schrodinger equation, and an 
alternative interpretation of the level-elimination at work in the 'supersymmetric' third 
spectral equivalence. 
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A Solving radial Schrodinger equations using Maple 

In this appendix we discuss the numerical treatment of the radial Schrodinger equation. 
The standard method is to integrate the ODE at varying values of energy, imposing 
boundary conditions either at the origin or infinity, and searching for the values of 
energy at which the other boundary condition is satisfied (see for example |^^). This 
approach runs into problems in the region ?R.el < —1/2, where the eigenvalue problem is 
better defined via analytic continuation. The series solution alternative that we describe 
here avoids this difficulty, can be implemented in only a few lines of Maple and seems to 
be a rather efficient method for finding the lower-lying levels, at least for a polynomial 
potential. We chose not to use a series of the sort described in §8 above, as the factor of 
exp(— means that any finite truncation of the series always decays at x — > oo; it is 
only in the infinite sum that the exponential growth of a solution at generic E is recov- 
ered. This makes it hard to detect eigenvalues reliably. Instead, we generated a pure 
power series directly in Maple, using an algorithm based on the method of Cheng ||3^ . 
The power series y produced by the program depends both on x and on E, and by con- 
struction it satisfies the boundary condition at x=0. At an eigenvalue, the solution must 
decay at large x, and by choosing a suitably-large value xq and searching for values of 
E at which y(xo,E)=0, the eigenvalues can be located with high accuracy. The value 
of Xq must be large enough that the asymptotic behaviour of the true solution has set 
in, but small enough that the approximated power series can be relied on. (The level 
of the approximation is controlled by the variable iterations in the program below.) 
This can be checked by examining plots of the candidate wavefunctions. We first give 
the code that we used, the particular example producing figure |l] below. The values of 
a and / are specified in the second line. 

> Digits: =20: iterations : =40 : 

> alpha: =-4: 1:=0: 

> V:=x~6+alpha*x~2; 

> L:=poly->suni(coeff (poly,x,n)*x~(n+2)/(n+2)/(n+2*l+3) ,n=0. .degree(poly,x)) : 

> P:=l:for i from 1 to iterations do P:= simplify (1+L( (V-E) *P) ) end do: 

> y:=x~(l+l)*P: 

> spectrum:=f solve(eval(y,x=3.2)=0) ; 

A specific level, determined by the integer plotlevel, is then plotted as follows: 

> with(plots) : plotlevel : =0 :Ee : =spectrum[plotlevel+l] ; 

> xmax:=2.5:ymin:=-35:ymajc:=100: 

> display([ 



20 



plot(eval(80*y,E=Ee) ,x=0. . xmax , ymin . .ymax,color=blue,linestyle=2,thickness=2) , 
plot ( (V+1* (l+l)/x~2) ,x=0. . xmax , ymin . .ymax,color=red,thickness=2) , 
seq(plot( [[0.05, spectrum [lev]] , [0. 12 , spectrum [lev] ] ] , 

color=black, linestyle=l ,thickness=3) , lev=l. .7), 
plot ([[0.01, spectrum [plotlevel+1] ] , [0. 17, spectrum [plotlevel+1] ] ] , 
color=black, linestyle=l ,thickness=l) 

]); 

The levels are contained in the list spectrum, and in table |^ below we compare semi- 
classical, nonlinear integral equation (NLIE) and Maple results for two sides of the 
5C/(3)-inspired duality of §5. To compile the table, we increased Digits to 40 and 
iterations to 50; nevertheless, each column of data still took less than 6 minutes of 
CPU time on a 650 MHz Pentium III machine, running under Linux. 





(a = 0.6,/ = 0.3) 


{a = 2.1, 


I = 0.05) 


Esc 




-S'Maple 


-£^NLIE 


Eyiaplc 


5.9321 


5.968835071930586 


5.968835071930586 


5.968835071930586 


5.968835071930586 


17.4454 


17.45222031394444 


17.45222031394444 


17.45222031394444 


17.45222031394444 


32.4611 


32.46443928842889 


32.46443928842889 


32.46443928842892 


32.46443928842889 


50.2951 


50.29711627959066 


50.29711627959066 


50.29711627959066 


50.29711627959066 


70.5566 


70.55798308714998 


70.55798308714998 


70.55798308714998 


70.55798308714998 


92.9834 


92.98438725677185 


92.98438725677185 


92.98438725677185 


92.98438725677184 


117.3836 


117.3842907593505 


117.3842907593502 


117.3842907593505 


117.3842907593506 


143.6087 


143.6092804003097 


143.6092804003098 


143.6092804003097 


143.6092804003097 


171.5397 


171.5401630800539 


171.5401630800552 


171.5401630800539 


171.5401630800544 


201.0781 


201.0784630892047 


201.0784630892253 


201.0784630892047 


201.0784630892116 


232.1408 


232.1410302822398 


232.1410302455113 


232.1410302822398 


232.1410302776264 


264.6562 


264.6564313215198 


264.6565101370171 


264.6564313215198 


264.6563953856096 


298.5623 


298.5624300607573 


299.0784472732795 


298.5624300607573 


298.7522851488036 


333.8040 


333.8041674124193 


339.1746343166020 


333.8041674124193 


347.4041682854545 



Table 1: Numerical Results 



From the table, it is clear that the NLIE ( p. 15 ) is able to find the energy levels with 
high accuracy, provided |a|<M-|-l — |2/-|-1|, and for such values of the parameters 
it seems to be the most reliable method to find the full set of energy levels. The power 
series approach also works extremely well for the low-lying energy levels, but loses 
accuracy for the higher levels, at least if the value of iterations is kept reasonably 
small. Nevertheless, combining power series and semiclassical methods allows us to 
obtain good results, over the full spectrum, for any values of a and /. 
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The figures which follow illustrate some typical wavefunctions. In each case the 
low- lying energy levels are shown as bars near the left-hand side of the plot, with the 
level under scrutiny having double length. The dashed line is the (un-normalised) 
wavefunction associated with this energy level, and the solid line the potential itself. 

The first set of four figures illustrates the second spectral equivalence, connecting the 
potentials with a, I = —4,0 and 7/2, —5/4 respectively. Note that the dual problem is 
an example where the eigenfunction is not square-integrable; nevertheless, the numerical 
solution converges and reproduces exactly the required spectrum. Figures |^ and |2| show 
the ground state, and figures ^ and ^ the first excited state. 




2.5 



Figure |: a = -4, / = 0, Eq = 1.0057683 Figure |: a = 7/2, I = -5/4, Eq = 1.0057683 




2.5 



Figure |: a = -4, / = 0, Ei = 10.572585 Figure |: a = 7/2, I = -5/4, Ei = 10.572585 
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Figure ^ depicts the quasi-exactly solvable case for the third energy level, at J = 2 
and / = 0. This corresponds to the ground state of the SUSY (H) dual potential, shown 
in figure ^ The fourth level of the QES problem has the same energy as the first excited 
state of the SUSY dual, and these two are shown in figures and Figures ^ and |l^ 
then show the first two states of the potential related to the SUSY dual by the second 
spectral equivalence, the tilde-duality T. 




2.5 



Figure |: a = -9, / = 0, E2 = 16.919850 Figure ^. a = 3, I = 2, Eq = 16.919850 




2.5 



Figure 0: a = -9, / 0, E3 = 32.240265 Figure ^. a = 3, I = 2, Ei ^ 32.240265 



23 





Figure |: a = 6, / = 3/2, Eq = 16.919850 



Figure 0: a = 6, / = 3/2, = 32.240265 



Finally, in figures 11 and |l^ we illustrate how 'extra' energy levels can appear with 
the irregular boundary condition in resonance situations, a phenomenon that was men- 
tioned at the end of the conclusions above. Again, we take J = 2. To avoid numer- 
ical difficulties, in figure |l^ we shifted the angular momentum slightly away from the 
exactly-resonant value. 




Figure O: Q = -9, / = 0, Eq = -4.898979 




Figure |1J: a = 6, I 
Eo = -4.898974 



5/2+10- 
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An elementary proof of the Bessis, Zinn-Justin, Bender 
and Boettcher conjecture 



B 



A conjecture of Bessis and Zinn-Justin |^2[, generalised by Bender and Boettcher ||2 
states that the eigenvalues of the PT-symmetric Schrodinger equation 



[IX 



,2M 



(B.l) 



are real and positive for M > 1. The contour C on which the wavefunction is defined 
can be taken to be the real axis for M < 2 ; beyond this point, the contour should 
be deformed down into the complex plane so as to remain in the same pair of Stokes 
sectors p3|. (For an informal review in the context of the ODE/IM correspondence, 
see also [0.) This conjecture has provoked a fair amount of work in recent years, a 
sample being refs. [24-|3^. In this appendix we consider a slightly more general class of 



"PT-symmetric spectral problems, namely 



dx'^ 



[IX 



,2M 



a{ix) 



M-l 



Kl + l) 



iJk{x) = Xk Mx) , Mx) G L\C) (B.2) 



with M, a and I real. Again, for M < 2 the contour C can be taken to be the real axis, 
though if 1(1+1) 7^ it should be distorted so as to pass below the origin. We shall 
prove reality of the spectrum for M > 1, a < M+l+|2^+l|, and positivity for M > 1, 
a < M+1— 12/+1|. The spectrum might be real for a greater range of a, but strict 
positivity certainly fails on the lines a = M+l— |2Z+1|. Even with the restrictions 
on a, our result includes the previously-considered cases: for a = 1(1+1) = and 
M = 3/2, a version of the original Bessis - Zinn-Justin conjecture is recovered; allowing 
M to vary then gives the generalisation discussed by Bender and Boettcher, while the 
conjecture for a = and I small was proposed in [0]. (Strictly speaking the original 
BZ-J conjecture concerned the potential x'^ + igx^ with g real; our discussion applies to 
the strong-coupling limit of this problem.) 
Setting <^(x) = ip(x/i), ( [B.2[ ) becomes 



dx^ 



+ x^^' + ax^''' + ^-^fi 



1>k(x) = -Xk<^k(x), <l>k(x) e L\iC) , (B.3) 



and has the same form as (2.1) with E = —Xk, though with different boundary condi- 
tions: to qualify as an eigenfunction, <I> must decay as — > oo along the contour iC. 



However, it is an easy generalisation of the discussion in §7 of [11 1 that the function 
T(—X,a,l) defined in (2^) is the spectral determinant associated to the spectral prob- 
lem ( |B.3| ). This identification allows us to study the generalised BZ-JBB conjecture 
( |B.2| ) using techniques inspired by the Bethe ansatz. 
We start from equation ( 2.10| ): 



T^+\E)D^+\E) = l^-(2'+l+-)/2D(-)(172Af^) + ^(2i+l+a)/2^(-)(^-2M^) ^ (g^4) 
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and define the zeroes of T^'^\E) = T(E, a, I) to be the set {—Xk}- (Note that for a = 0, 
(B.4) reduces to the T-Q system obtained in Putting E = —X^ in ( |B.4 ) and using, 



for M > 1, the factorised form for \E) gives the following constraints on the X^s: 



En ^ + ^ "^^^ Xk 



n I 1=-^''''^'"^ k=o,i,... (B.5) 



n=0 



Since the original eigenproblem ( [B.2| ) is invariant under / — > —l—l, we can assume 
/ > —1/2 without any loss of generality. Then each En ^ is an eigenvalue of an Hermitian 



operator 7i{M, —a, I), and hence is real. Furthermore a Langer transformation (see 
also [||,^) shows that the En ^ solve a generalised eigenproblem with an everywhere- 
positive 'potential', and so are all positive, for a < 1+21. This can be sharpened by 
considering the value of D^~\E)\^^^^. From ( [B.15| ) below, this first vanishes when 

a = M+21+2. Until this point is reached, no eigenvalue En ^ can have passed the 
origin, and all must be positive. (It might be worried that negative eigenvalues could 
appear from E = —oo , but this possibility can be ruled out by a consideration of the 
Langer-transformed version of the equation.) 

Taking the modulus^ of (jB.Sp , using the reality of the E^ , and writing the eigen- 



values of ( [B.2D as Afc = |Afc| exp(i5fc), we have 

- nEi^)^ + \Xk\' + 2Ei-^\Xk\cos{j^,+6k) \ ^ 
[{Ei^y + lAfeP + 2Ei-^\Xk\ cos(^ -6k)) ' 

For a < M + 2Z + 2 , all the En ^ are positive, and each single term in the product on 



the LHS of ( |B.6| ) is either greater than, smaller than, or equal to one depending only 
on the relative values of the cosine terms in the numerator and denominator. These are 
independent of the index n. Therefore the only possibility to match the RHS is for each 
term in the product to be individually equal to one, which for / requires 

cos{j^ + Sk) = cos{j^ - 6k) , or sin( sin(5fc) = . (B.7) 

Since M > 1, this latter condition implies 

5k = mr , n G Z (B.8) 



and this establishes the reality of the eigenvalues of (|B.2| ) for M > 1 and a < M + 21 + 2 
or, relaxing the condition on /, a < M + 1 + |2^-|-1|. 

One might ask what goes wrong for M < 1, since from (and, for the case 
I 7^ 0, |]lT]) it is known that most of the become complex as M falls below 1, at 
least for a = 0. The answer is that if M < 1, the order of D^~\E) is greater than 1, 
the factorised form of D^~\E) provided by Hadamard's theorem no longer has such a 
simple form, and the proof breaks down. 
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The borderline case M = 1 is the simple harmonic oscillator, exactly solvable for all 



/ and a . Starting from the discussion in §3 of [11|, it is easily seen that 



27r 



2 ' 4 / 4 

and so the eigenvalues of ( [B.2| ) are at A = 4n + 2 — ait (2/+1), n = 0, 1, . . . . All are 
real for all real values of a and and all are positive for a < 2 — |2/+1|. 

To discuss positivity at general values of M > 1, we can continue in M, a and I 
away from a point in this latter region, {Af=l,a < 2 — |2/+1|}. So long as a remains 
less than M + 1 + |2Z+1|, all eigenvalues will be confined to the real axis during this 
process, and the first passage of an eigenvalue from positive to negative values will be 
signalled by the presence of a zero in r(— A, a, at A = 0. Fortunately, T(— A, a, OIa=o 
can be calculated exactly, extending an argument given for a = in First, one 

notices that the function 

M+Q M — 1 / 1 2 \ 

ip{x) = ( 2M+2 y / (4^hi)M+i xM+T ,^,a,;J (B.IO) 



solves the Schrodinger equation 

, 2 _^2^ , 7(7 + 1) 



+ X — (TX M+1 + 



(^(j;) = A , (B.ll) 



where 

2M 



2Z+1 1 ^_ 2a 
M+1 ~ 2 ' ~ ~M+1 

(This transformation, which can be found via a pair of Langer transformations, leads 
to equation (pT^) in the case M=3.) Further, (/?(x) has the large-x asymptotic 



(fix) ~ ^ exp ) . (B.13) 



At i?=0, (7 = and ( B.llD is the simple harmonic oscillator, which can be solved exactly 



in terms of the confluent hypergeometric function U{a,b,z). Matching asymptotics at 
large x, 

^{x)\^^, = 4= e--'/' U ( i(7+i)-iA, 7+i, ) . (B.14) 

Reversing the variable changes, extracting the leading behaviour as a; ^ and compar- 
ing with (p. 9]) , we find 

^ 2i+i-c 1 r(^] 

D(E,aJn,__„ = = — J^., . (B,15) 

^ I 2M+2 2 
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Now T{E,a,l)\j^^Q follows from p^ ), remembering that D^-\E) = D{E,-a,l) : 
T{E,a,l)\j,_, = T(+\E)\j,_,= (M±l)^' ^ ^ . (B.16) 



^ 2 ^ 2M+2 J \2 2M+2 

The first zero arrives at = —A = when a = M + 1 — \2l+l\, and so for all 
a < M + 1 — |2Z+1| , the spectrum is entirely positive, as claimed. 

In finishing, we return to the reality of the spectrum encoded by T(— A). We have 
proved that, if M > 1, the eigenvalues Afc are real for all real a < M + 1 + |2/+1|. 
One might conjecture that this reality should hold for all real a and I. However, this is 
definitely not the case: for M = 3, at the QES points and with / sufficiently negative, an 
examination of the Bender-Dunne polynomials shows that the exactly-calculable part 
of the spectrum of D{E,aj) has at least one pair of complex-conjugate eigenvalues. 
The identity T{—E,Maj) = 'j{aj)D{E, aj), which follows from the results obtained in 
§7 above, shows that T(—E,Maj) must share these complex zeroes. If {aj,l) are real 
then so are (Sj, Ij), and so such examples demonstrate that T can have complex zeroes 
even while M > 1 and a and / are real. It would be worthwhile to map out the full 
extent of the region where the spectrum is entirely real, but we have not yet done this, 
beyond a quick check that it appears to extend at least some way beyond the domain 
a < M + 1 + |2Z+1| covered by the proof given in this appendix. 

Note added in proof. 

(1) We have recently obtained some further results on the region within which the spectrum 
of the T^T-symmetric problem discussed in appendix B becomes complex. These can be found 
in @. 

(2) An alternative treatment of a class of T^T-symmctric quantum mechanical problems similar 
to those discussed in appendix B can be fomid in H, Ea . 
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